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Abstract — This paper presents a novel technique for 
designing an Infinite Impulse Response (IIR) Filter with 
Linear Phase Response. The design of IIR filter is always a 
challenging task due to the reason that a Linear Phase 
Response is not realizable in this kind. The conventional 
techniques involve large number of samples and higher 
order filter for better approximation resulting in complex 
hardware for implementing the same. In addition, an 
extensive computational resource for obtaining the inverse 
of huge matrices is required. However, we propose a 
technique, which uses the frequency domain sampling along 
with the linear programming concept to achieve a filter 
design, which gives a best approximation for the linear 
phase response. The proposed method can give the closest 
response with less number of samples (only 10) and is 
computationally simple. We have presented the filter design 
along with its formulation and solving methodology. 
Numerical results are used to substantiate the efficiency of 
the proposed method. 

Index Terms — Infinite impulse response, linear 
programming, filter design, linear phase filters 

I. Introduction 

The design of the infinite impulse response filter is 
often restricted to the specification of the magnitude 
response due to the fact that a linear phase response is not 
physically realizable. Some methods in the literature try 
to achieve a good approximation of the linear phase 
response, but results in extensive and complex hardware. 
Thus, the design phase involves huge computational 
resources as well as implementation phase requires 
extensive hardware. We intend to simplify the 
computation in design as well as minimize the hardware 
by just opting for a second order filter. We introduce a 
concept of using the frequency domain samples of the 
desired magnitude and phase response of the filter and 
using the same to obtain the filter coefficients to 
implement the filter. The digital filter can be modeled as 
a linear system represented by 



Ax = b 



(i) 



Where, A is the matrix containing the sampled frequency 
responses (magnitude & phase) of the given filter [1]. 
The given magnitude and phase responses of the filter 
can be sampled at a convenient sampling frequency and 
the number of samples can be conveniently chosen. We 
will now formulate the filter design problem as a linear 
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programming problem and then we will redefine the same 
problem with some constraints for solving the linear 
programming model. The frequency domain responses of 
the filter (both magnitude and phase) are used to model 
the linear programming problem. Let us consider the 
standard expression [2] for the filter transfer function, 



Y(Z) b +b lZ - 



U(z) 



l + a_z~ 



-...a n z 



_v < 

-D 



(2) 



Where U (z) is the z-transform of the input signal, Y (z) is 
the z-transform of the output signal, and 'a' and 'b' 
factors are real valued coefficients. From the above 
transfer function, we can write the time-domain 
difference equation that implements the filter as, 



y(k) = -a l y(k-l)-....a D y(k 

+ b u(k) + ... + b N u(k - N). 



D) 



(3) 



where the a and b coefficients are exactly the same as in 
the transfer function above, k is the time index, u (k) and 
y (k) are the current values of the input and output 
(respectively), u(k- N) was the input value N samples in 
the past, and y(k- D) was the output value D samples in 
the past. Generally, if the frequency response of a system 
at a frequency U) x is given in magnitude/phase form 
as A l /-(f) l the output amplitude will be A x times the input 
amplitude and the output phase will be shifted an angle 
(pi relative to the input phase when a steady-state sine or 
cosine wave of frequency cul is applied to the system. 
For instance, if the input to the system described above at 
time instance 'k' is u x (k) = Cos(kiOits), where t s is the 
sampling period, then the output will be 
y t (k) = A, Cos(ku>i ts + q>\). The input and output values 
at any sample time can be determined in a similar 
manner. The final sampled values in the matrix form is, 



u,(-N) 
u 2 (-N) 



M0) 




-y,(-i) •• 

-_(-!) .. 


- _(-D) 


«,(0) 
« 2 (0) 


yAO) 




->v(-i) ■• 


-y u (-D) 


M0) 



u u (-N) 

(4) 

The above equation is similar to the system A x = b ; as 
given by (1). In the above formulation, we have set the 
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current sample index as zero. Our aim now is to solve the 
above problem, and get the filter coefficients given by the 
vector 'x'. Standard methods use the inverse of the 
matrix 'A' which involves a heavy computation for huge 
matrices [3], [4], [5], [6], [7], [8]. If the inverse of the matrix 
does not exist, then Singular Value Decomposition can be 
used for Pseudo inverses [9]. To avoid this and make the 
procedure computationally simple, we introduce the 
concept of linear programming into this context. We now 
formulate the above as a simple optimization problem 
with few constraints. 

II. Problem Formulation 

In the above model discussed, the matrix 'A' contains 
the sampled coefficients at different frequencies and the 
column vector 'x' contains the filter coefficients. The 
first 'D' values in 'x' represent the coefficients ai, a 2 , etc 
in denominator and the next 'N' values give the b , bi, b 2 , 
etc coefficients in numerator. The column vector 'y' is 
the desired output value (both amplitude and phase) at 
specific frequencies. Now we can move on to formulate 
the optimization problem to find the filter coefficients. 
We first aim to find the best possible solution by 
introducing a factor 's', which is vector containing 
allowable error. Then we try to reduce the error to get the 
best solution by minimizing this error parameter. 



Ax + £ = y; 
=^> £ = y - Ax 



(5) 



Now our goal is to frame a minimization problem to 
reduce the error. Here our variables are 'x' and '£ '. 



Min £>,+ £ 



a : 



S.t Ax + £ = y; < £ < a.*y anda>0; 



(6) 



Instead of directly manipulating the error parameter, we 
can define two variables Wi and w 2 (essentially vectors of 
length same as 'e') such that 



£ = Wj -w 2 ; 

So our equation (5) becomes, 

Ax + w t -w 2 =y; 



(7) 



(8) 



We also have a possibility to restrict the maximum value 
error to a limit say a fraction of the actual response value, 
so that relatively there is no adverse change in the overall 
response. Letting that positive factor to be 'a' that 
contains the fraction (eg: 0.01 or 1/100 th) , where 'a' is 
also a vector. 



0<£<a.*y 



(9) 



Here '.*' is a Matlab operation, which is an element by 
element multiplication. Rewriting the above constraint 
using Wj and w 2 , 

0< Wj <a.*y 

0<w 2 <a.*y (io) 

Now the problem, which we need to solve, is given by 



Min — 2J* X +— Z, w 2+Z, a i 

X,W x ,W 2 ,a <-l i ^2 i i 

s.t Ax + Wj-w 2 =y; 

0<Wj<a.*y, 0<W 2 <a.*y anda>0; 

(11) 

Here Ci and C 2 are control parameters that can be 
adjusted to vary the depth of optimization of Wi and w 2 
respectively. By the above formulation, we try to reduce 
the value of error 'g' along with the maximum 
permissible limit defined by 'a'. When we attempt to 
reduce the value of 'a', indirectly the maximum value 
that Wi and w 2 can assume is reduced which in turn 
results in the value of '£'. The problem is a chain reaction 
one. That mutually reduces the above said parameters. 
This formulation of the optimization problem is easily 
solvable by any simple tool, even with a spreadsheet 
package like Microsoft Excel. This formulation results in 
reduced computational complexity and less number of 
frequency domain samples (eg: 8 or 16). 

III. Numerical Results 

We have taken few examples [10], [11], [12] of low 
pass; band pass and high pass filter and attempted to 
solve for the filter coefficients. The knowledge of desired 
frequency response (both magnitude as well as phase) is 
the only requirement. We have considered a second order 
filter in all the cases discussed with 10 sample values. 
The number of samples is left to the choice of the 
designer with only condition that minimum number of 
samples be 10 for a good approximation. The proper 
choices of sample points play an important role in the 
response approximation. The sampling is done at the rate 
of 1000Hz or 10" 3 seconds. Since we have chosen a 
second order filter, we have N=D=2. We have solved 
the optimization problem using Microsoft Excel and the 
results are verified by using Matlab. The obtained 
coefficients are used to plot the magnitude and phase 
response and the results are verified. 

A. Illustration 1: 

A Low pass IIR Filter was designed for the given 
Specifications as shown in figure 1. The signal is 
sampled and the sample values corresponding to 
amplitude and phase are tabulated in Table 1. The filter 
problem is formulated as a linear problem using the 
above discussed formulation and is solved using Excel 
and the coefficients obtained are listed in Table 2, also 
the response obtained using Matlab is shown in figure 1. 
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TABLE 1: SAMPLED DATA :( 10 POINTS) - LPF 



Amplitude 


Phase (radians) 


Frequency (Hz) 


1.0 


-0.1 


2.7400 


1.0 


-0.5 


50.0400 


1.0 


-0.97 


97.9600 


1.0 


-1.01 


101.8700 


0.9 


-1.1 


110.0000 


0.6 


-1.25 


125.0000 


0.5 


-1.25 


150.0000 


0.3 


-1.25 


178.2800 


0.1 


-1.25 


300.4200 


0.001 


-1.25 


450.5600 



TABLE 2: OBTAINED COEFFICIENTS - LPF 



a 


ai 


a 2 


bi 


b 2 


b 3 


1.00 


-0.8663 


0.3141 


0.2141 


0.1122 


0.1209 



It can be inferred that the phase response exhibits a linear 
characteristics in our pass band which is - 100 Hz. The 
phase is exactly '-1 radians' at 100 Hz and it exactly fits a 
straight line with slope angle 45 degrees. We are not 
concerned of the phase relation in the stop band, but the 
phase seems to fit a linear curve in stop band also. The 
small deviation in the response in the transition from the 
pass band to stop band can be rectified by increasing the 
number of sample points in this region. We have 
experimented with 10 samples and tried to prove that 
with just 10 samples, we are able to get a good 
approximation. Therefore, by suitably increasing the 
number form 10 to 16 or 32, the obtained response will 
exactly follow the desired one. 




D 
-03 




T .'v-.i Response 










/' 


De.ircd 

Obtained 








04 

■0.6 
J™ 


/ 






t 

t 


i 










7 




1 -12 






/ 


a. 

-1 4 

-1.6 

-i.a 




t 
/ 


• 













50 100 150M02S>3tB350MO.I£i]H>0 

Frequency 

Fig lb. Comparison of desired and obtained Phase Responses of the low 
pass filter 

The desired and obtained magnitude and phase responses 
are comparatively studied and the figure 1 depicts the 
same. It is quite clear that the obtained response better 
approximates the desired one with just 10 samples and a 
second order IIR filter. With this method, a linear phase 
response can be achieved even with an IIR filter. 

B. Illustration 2: 

We now discuss the second example of designing a 
band pass filter to operate in the band from 100 Hz 200 
Hz with a linear phase response as depicted in fig 2b. The 
signal is sampled and the sample values corresponding to 
amplitude and phase are tabulated in Table 3. The filter 
problem is formulated as a linear problem using the 
above discussed formulation and is solved using Excel 
and the coefficients obtained are listed in Table 4 and the 
response obtained using Matlab is shown in figure 2. 

TABLE 3: SAMPLED DATA (10 POINTS) - BPF 



Amplitude 


Phase (radians) 


Frequency 
(Hz) 


0.0 


0.5 


0.0000 


0.25 


0.5 


50.0000 


0.5 


0.5 


95.0000 


0.9 


0.5 


100.0000 


1.0 


-0.11 


150.0000 


1.0 


-0.63 


200.0000 


0.75 


-1.0 


250.0000 


0.5 


-1.46 


300.0000 


0.3 


-1.5 


350.0000 


0.2 


-1.5 


450.5600 



Fig la. Comparison of desired and obtained Magnitude Responses of the 
low pass filter 
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TABLE 4: OBTAINED COEFFICIENTS - BPF 



220 

a 


ai 


a 2 


bi 


b 2 


b 3 


1.00 


-0.7502 


0.2644 


0.2617 


0.2471 


-0.5010 



The comparison plots (figure 2) between the desired and 
obtained responses are presented. The phase response in 
the pass band, that is 100 Hz - 200 Hz, is quite linear 
with respect to frequency. This region is of interest to us 
and is the key feature to be noticed. 



Magnitude Response 




Fig 2a. Comparison of desired and obtained Magnitude Responses of 
the Band pass filter 



Phase Response 




Fig 2b. Comparison of desired and obtained Phase Responses of the 
Band pass filter 

C. Illustration 3: 

The Third example we have considered is that of a 
high pass filter. The procedure similar to the previous 
examples is followed here also. The Desired responses 
(figure 3), sampled data (table 5), the obtained 
coefficients (table 6) and obtained responses (figure 3) 
are presented below. The comparison plot (figure 3) 
clearly describes the linear phase characteristics of the 
designed filter. By increasing the number of samples, we 
can make the obtained response exactly follow the 
designed one. Here we have focused on the linear phase 
response, hence chosen just 10 samples. By choosing 32 
samples, the obtained response is found to be the exact 
copy of the desired response. 



TABLE 5 


SAMPLED DATA (10 


JOINTS) - HPF 


Amplitude 


Phase (radians) 


Frequency (Hz) 


0.0 


1.7 


0.0000 


0.1 


1.7 


50.0000 


0.2 


1.7 


95.0000 


0.3 


1.7 


100.0000 


0.5 


1.5 


150.0000 


1.0 


1.2 


200.0000 


1.0 


0.9 


250.0000 


1.0 


0.6 


300.0000 


1.0 


0.4 


350.0000 


1.0 


0.1 


450.5600 



TABLE 6: OBTAINED COEFFICIENTS - HPF 



a 


ai 


a 2 


bi 


b 2 


b 3 


1.00 


-0.2194 


0.1662 


0.4975 


-0.7182 


0.2208 



Magnitud* Respante 




Fig 3a. Comparison of desired and obtained Magnitude Responses of 
the high pass filter 



Phase Response 




Fig 3b. Comparison of desired and obtained Phase Responses of the 
high pass filter 

IV. Parameters Influencing The Filter Design 

Firstly, the filters designed by the above method solely 
depend upon the number of samples and the sample 
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values. The key point is the choice of the sample values. 
The frequency points where there is a significant change 
in the response either magnitude or the phase should be 
captured by the sample values we choose. For a normal 
filter, we have considered (low pass or high pass or band 
pass); four points in the transition band, three points in 
the pass and stop band will do the trick. However, if the 
filter should provide specific responses, then the number 
of samples should be chosen accordingly depending upon 
the level of approximation we require. However, as a 
general case a filter designed by this proposed 
methodology with 10 samples should yield a good 
approximation of the magnitude and phase response 
when compared to the design by existing methods. In 
addition, we can achieve a linear phase Infinite Impulse 
Response (IIR) filter which is normally not realizable. 
Hence, this method can be effectively used to achieve the 
desired response using an IIR filter. 

Secondly, we are interested only in the phase response 
in the pass band and the phase in the stop band is not of 
our concern. Yet, the literatures reveal that the phase 
change in stop band will influence the overall response of 
the filter response. These gives a cushion of changing the 
phase in the stop band and further achieve better results. 
Therefore, we can introduce few more variables that are 
phase values in stop band in our formulation, even 
provide some lower and upper bounds, and further refine 
the formulation. 
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V. Conclusion And Future Scope 

The formulation of the IIR filter design problem, with 
both the magnitude and phase specifications, is easily 
done using this proposed technique of linear 
programming. This application of linear programming for 
design of linear phase IIR filters is a novel methodology. 
This requires only very low computational resources in 
spite of a very good approximation of the linear phase 
response. The results obtained by the proposed method 
are quite encouraging. The results can be further 
improved by proper choice of the number of samples and 
position of sample points. These results presented above 
guarantee the efficiency of the proposed method. The 
literature shows that by slightly tuning the phase in the 
stop band (which is of no significance) the phase 
response in the pass band can be changed accordingly. 

These phenomena can be incorporated in the proposed 
technique to achieve best results. The combination is 
implemented and the performance is being investigated. 
It is very much conducive that the proposed method may 
be effectively adopted by the design engineers. The 
proposed method proves to be a simple method when 
compared to the counterparts for providing better results. 
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